Leptonic decay constants fr> and /d in three flavor lattice QCD 
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We determine the leptonic decay constants fn, and fo in three flavor unquenched lattice QCD. We use 
0(a 2 )-improved staggered light quarks and 0(a)-improved charm quarks in the Fermilab heavy quark formalism. 
Our preliminary results, based upon an analysis at a single lattice spacing, are fo B = 263lg ± 24 MeV and 
f D = 225^13 ± 21 MeV - In each 

case, the first reported error is statistical while the second is the combined 

systematic uncertainty. 



1. INTRODUCTION 

The leptonic decay constants /b, and fs are 
critical in testing the flavor sector of the Standard 
Model. Reliable lattice calculations are of fun- 
damental importance since the determination of 
these decay constants remains beyond the reach 
of experiment. 

Precise experimental determinations of the de- 
cay constants fu 3 and fjj and the semileptonic 
decays D — > tt£i/ and D — > Klv will result from 
the high-statistics charm program of CLEO-c. 
Comparing these experimental results with lat- 
tice results will serve both as a critical check of 
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lattice methods for charm and as a means of as- 
sessing the reliability of lattice calculations for the 
bottom quark. 

This work calculates the leptonic decay con- 
stants fo s and fjj using 0(a 2 )-improved stag- 
gered light quarks and 0(a)-improved charm 
quarks in the Fermilab heavy quark formalism 
[1]. It was done with three flavors of light sea 
quarks. The staggered fermion action for the 
light quarks makes possible calculations at lighter 
quarks masses than previously used [2], allowing 
a better controlled chiral extrapolation. We use 
the results of staggered, partially quenched chi- 
ral perturbation theory (SxPT) in performing the 
chiral extrapolations. 
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2. f Dy /ME/fn e y/rnE; RATIO 

We use the MILC collaboration Asqtad gauge 
ensembles with lattice spacing a w 1/8 fm [3]. 
Each ensemble has one sea quark flavor of mass, 
m s , approximating the strange quark and two 
flavors, degenerate in mass, m u . The five en- 
sembles we use have light flavors in the range 
0.1m s < m u < 0.6m s . For each ensemble, de- 
cay constants were computed at twelve logarith- 
mically spaced values of the valence quark mass 
in the range 0.1m s < m q < m s . In all, decay con- 
stants were computed for sixty partially quenched 
(m q ,m u ) combinations. 

The chiral expansion for the decay constants is 
known at leading order in both the heavy quark 
expansion and staggered chiral perturbation the- 
ory [4]. We take the parameterization 

Rq/s = 1 + ao(A/ g - A/ a ) + (m, - m s ) x 

(ai + a 2 fh + a 3 m q + a 4 m s + . . .) (1) 

for the ratio R q/s = fD qy /fn^/f Dsy /fn^ in our 
chiral extrapolation. For staggered quarks, the 
chiral log terms Af x contain discretization effects 
from taste violations in the pseudoscalar masses 
as well as explicit taste violation terms. We in- 
clude quark mass terms up to 0(m 2 ) with the 
constraint 04 = a 3 — a 2 . The sum of sea quark 
masses is rh = 2m u + m s . 

The parameters aj in Eq. (1) are determined in 
a fully covariant x 2 minimization using all sixty 
partially quenched decay constant results. Pseu- 
doscalar masses and coefficients of the explicit 
taste violation terms were fixed to the values de- 
termined in an analysis of the light mesons [5]. 
Bayesian priors were input for each a,j. With 
£ = 1/(4tt/) 2 , the prior for a a is -0.5£(1 + 
3# 2 )(1 ± 0.30) with g 2 w 0.35. Other priors are 
± 1 in units of 2£s where s is the slope relating 
the quark mass and the pion mass squared. 

The fit, shown in Fig. 1, has xV^of = 0.2. 
The extrapolation according to S%PT is shown 
along the "full QCD" direction where the valence 
quark mass equals the sea quark mass. The pa- 
rameters in this extrapolation are determined us- 
ing all sixty partially-quenched points. We obtain 
R d/s = 0.833±°;^ 2 in the chiral limit. The error 
is the combined uncertainty from statistics and 
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Figure 1. The chiral extrapolation of the ratio 
fny/nT-D I ' fD B ^TTbD s - The ratio is determined by 
the solid curve where staggering effects have been 
removed. The dashed curve includes these dis- 
cretization effects. 



parameter prior estimates for the extrapolation 
function. Systematic effects from matching the 
lattice theory to QCD, the lattice spacing and the 
tuning of charm quark mass mostly cancel in the 
ratio. Residual discretization effects from light 
quarks, not removed by the extrapolation proce- 
dure, are estimated to be 4%. Our estimate is 
based on the size of known taste-breaking effects 
shown in Fig. 1 and is similar to taste-breaking 
effects found in fx and f v [5] . The ratio was de- 
termined using the nominal strange quark mass 
rather than the tuned value which leads to an 
uncertainty of 2%. 

3. fo. yftnjDl DETERMINATION 

We obtain /p.y'm^ by first interpolating to 
the tuned valence m s value obtained from the 
light mesons on the same gauge ensembles. Then, 
a mild sea quark extrapolation (sec Fig. 2) is 
needed to obtain fr> s . We extrapolate linearly 
to the tuned rh = (m u + rrid) /2 value obtained 
from the light mesons [5]. The combined statis- 
tical and extrapolation error is 3.3%. The uncer- 
tainty from the tuning of m s and rh is 1%. The 
tuning of the charm mass leads to a 4% error. 
The lattice spacing uncertainty is 2% [6]. The 
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Table 1 

Error budget as percentage of each quantity. 
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Figure 2. Extrapolation in the light sea quark 
mass for fD B ^mD B - The curves show a linear 
fit (soild) and the 68% confidence level statistical 
error bounds (dotted). 



dominant systematic uncertainty 7%, is from the 
mismatch between the lattice theory and QCD, 
as discussed in Ref. [7]. Our final results will in- 
clude an improved estimate of this uncertainty 
incorporating results from finer and coarser lat- 
tice spacings, which are now in progress. 

4. RESULTS 

Statistical and systematic uncertainties are 
summarized in Table 1. Our estimates of heavy 
quark matching effects and light quark discretiza- 
tion effects are based on results from a single lat- 
tice spacing. We will refine our error estimates 
and update our results once decay constants from 
additional lattice spacings are known. The heavy 
quark matching uncertainty can be reduced by in- 
cluding the higher order matchings for the action 
and the currents [8,9]. 

Combining in quadrature the systematic uncer- 
tainties shown in Table 1, we find our preliminary 
results: 




1.20 ±0.06 ±0.06 



263+^ ± 24 MeV , 
225±i3 ± 21 MeV 
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